Unprecedented insights into extents of biological responses to physical forcing in an Arctic sub-mesoscale filament by combining high-resolution measurement approaches.

In Fram Strait, we combined underway-sampling using the remote-controlled Automated Filtration System for Marine Microbes (AUTOFIM) with CTD-sampling for eDNA analyses, and with high-resolution optical measurements in an unprecedented approach to determine variability in plankton composition in response to physical forcing in a sub-mesoscale filament. We determined plankton composition and biomass near the surface with a horizontal resolution of ~ 2 km, and addressed vertical variability at five selected sites. Inside and near the filament, plankton composition was tightly linked to the hydrological dynamics related to the presence of sea ice. The comprehensive data set indicates that sea-ice melt related stratification near the surface inside the sub-mesoscale filament resulted in increased sequence abundances of sea ice-associated diatoms and zooplankton near the surface. In analogy to the physical data set, the underway eDNA data, complemented with highly sampled phytoplankton pigment data suggest a corridor of 7 km along the filament with enhanced photosynthetic biomass and sequence abundances of sea-ice associated plankton. Thus, based on our data we extrapolated an area of 350 km2 in Fram Strait with enhanced plankton abundances, possibly leading to enhanced POC export in an area that is around a magnitude larger than the visible streak of sea-ice.

Global Distribution and Diversity of Marine Parmales.

Parmales (Bolidophyceae) is a minor eukaryotic phytoplankton group, sister to diatoms, which exists as two distinct forms of unicellular organisms: silicified cells and naked flagellates. Since their discovery, many field studies on Parmales have been performed; however, their global distribution has not yet been examined in detail. We herein compiled more than 3,000 marine DNA metabarcoding datasets targeting the V4 region of the 18S rRNA gene from the EukBank database. By linking this large dataset with the latest morphological and genetic information, we provide updated estimates on the diversity and distribution of Parmales in the global ocean at a fine taxonomic resolution. Parmalean amplicon sequence variants (ASVs) were detected in nearly 90% of the samples analyzed. However, the relative abundance of parmaleans in the eukaryotic community was less than 0.2% on average, and the estimated true richness of parmalean ASVs was approximately 316 ASVs, confirming their low abundance and diversity. A phylogenetic ana-lysis divided these algae into four clades, and three known morphotypes of silicified cells were classified into three different clades. The abundance of Parmales is generally high in the poles and decreases towards the tropics, and individual clades/subclades show further distinctions in their distribution. Collectively, the present results suggest clade/subclade-specific adaptation to different ecological niches.

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing.

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

Transcriptome analysis of the harmful alga Heterosigma akashiwo under a 24-hour light-dark cycle.

The photoperiod, which is defined as the period of time within a 24-hour time frame that light is available, is an important environmental regulator of several physiological processes in phytoplankton, including harmful bloom-forming phytoplankton. The ichthyotoxic raphidophyte Heterosigma akashiwo is a globally distributed bloom-forming phytoplankton. Despite extensive studies on the ecological impact of H. akashiwo, the influence of the photoperiod on crucial biological processes of this species remains unclear. In this study, gene expression in H. akashiwo was analyzed over a 24-hour light-dark (14:10) treatment period. Approximately 36 % of unigenes in H. akashiwo were differentially expressed during this 24-hour treatment period, which is indicative of their involvement in the response to light-dark variation. Notably, the number of differentially expressed genes exhibited an initial increase followed by a subsequent decrease as the sampling time progressed (T0 vs. other time points). Unigenes associated with photosynthesis and photoprotection reached their peak expression levels after 2-4 h of illumination (T12-T14). In contrast, the expression of unigenes associated with DNA replication peaked at the starting point of the dark period (T0). Furthermore, although several unigenes annotated to photoreceptors displayed potential diel periodicity, genes from various photoreceptor families (such as phytochrome and cryptochrome) showed unique expression patterns. Collectively, our findings offer novel perspectives on the response of H. akashiwo to the light-dark cycle, serving as a valuable resource for investigating the physiology and ecology of this species.

A single Prochlorococcus ecotype dominates the tropical Bay of Bengal with ultradian growth.

The Bay of Bengal (BoB) spans >2.2 million km2 in the northeastern Indian Ocean and is bordered by dense populations that depend upon its resources. Over recent decades, a shift from larger phytoplankton to picoplankton has been reported, yet the abundance, activity, and composition of primary producer communities are not well-characterized. We analysed the BoB regions during the summer monsoon. Prochlorococcus ranged up to 3.14 × 105 cells mL-1 in the surface mixed layer, averaging 1.74 ± 0.46 × 105 in the upper 10 m and consistently higher than Synechococcus and eukaryotic phytoplankton. V1-V2 rRNA gene amplicon analyses showed the High Light II (HLII) ecotype formed 98 ± 1% of Prochlorococcus amplicons in surface waters, comprising six oligotypes, with the dominant oligotype accounting for 65 ± 4% of HLII. Diel sampling of a coherent water mass demonstrated evening onset of cell division and rapid Prochlorococcus growth between 1.5 and 3.1 div day-1, based on cell cycle analysis, as confirmed by abundance-based estimates of 2.1 div day-1. Accumulation of Prochlorococcus produced by ultradian growth was restricted by high loss rates. Alongside prior Arabian Sea and tropical Atlantic rates, our results indicate Prochlorococcus growth rates should be reevaluated with greater attention to latitudinal zones and influences on contributions to global primary production.

Taxonomic difference in marine bloom-forming phytoplanktonic species affects the dynamics of both bloom-responding prokaryotes and prokaryotic viruses.

The production of dissolved organic matter during phytoplankton blooms and consumption by heterotrophic prokaryotes promote marine carbon biogeochemical cycling. Although prokaryotic viruses presumably affect this process, their dynamics during blooms are not fully understood. Here, we investigated the effects of taxonomic difference in bloom-forming phytoplankton on prokaryotes and their viruses. We analyzed the dynamics of coastal prokaryotic communities and viruses under the addition of dissolved intracellular fractions from taxonomically distinct phytoplankton, the diatom Chaetoceros sp. (CIF) and the raphidophycean alga Heterosigma akashiwo (HIF), using microcosm experiments. Ribosomal RNA gene amplicon and viral metagenomic analyses revealed that particular prokaryotes and prokaryotic viruses specifically increased in either CIF or HIF, indicating that taxonomic difference in bloom-forming phytoplankton promotes distinct dynamics of not only the prokaryotic community but also prokaryotic viruses. Furthermore, combining our microcosm experiments with publicly available environmental data mining, we identified both known and novel possible host-virus pairs. In particular, the growth of prokaryotes associating with phytoplanktonic organic matter, such as Bacteroidetes (Polaribacter and NS9 marine group), Vibrio spp., and Rhodobacteriales (Nereida and Planktomarina), was accompanied by an increase in viruses predicted to infect Bacteroidetes, Vibrio, and Rhodobacteriales, respectively. Collectively, our findings suggest that changes in bloom-forming species can be followed by an increase in a specific group of prokaryotes and their viruses and that elucidating these tripartite relationships among specific phytoplankton, prokaryotes, and prokaryotic viruses improves our understanding of coastal biogeochemical cycling in blooms.IMPORTANCEThe primary production during marine phytoplankton bloom and the consumption of the produced organic matter by heterotrophic prokaryotes significantly contribute to coastal biogeochemical cycles. While the activities of those heterotrophic prokaryotes are presumably affected by viral infection, the dynamics of their viruses during blooms are not fully understood. In this study, we experimentally demonstrated that intracellular fractions of taxonomically distinct bloom-forming phytoplankton species, the diatom Chaetoceros sp. and the raphidophycean alga Heterosigma akashiwo, promoted the growth of taxonomically different prokaryotes and prokaryotic viruses. Based on their dynamics and predicted hosts of those viruses, we succeeded in detecting already-known and novel possible host-virus pairs associating with either phytoplankton species. Altogether, we propose that the succession of bloom-forming phytoplankton would change the composition of the abundant prokaryotes, resulting in an increase in their viruses. These changes in viral composition, depending on bloom-forming species, would alter the dynamics and metabolism of prokaryotes, affecting biogeochemical cycling in blooms.

Phenology and ecological role of aerobic anoxygenic phototrophs in freshwaters.

BACKGROUND: Aerobic anoxygenic phototrophic (AAP) bacteria are heterotrophic bacteria that supply their metabolism with light energy harvested by bacteriochlorophyll-a-containing reaction centers. Despite their substantial contribution to bacterial biomass, microbial food webs, and carbon cycle, their phenology in freshwater lakes remains unknown. Hence, we investigated seasonal variations of AAP abundance and community composition biweekly across 3 years in a temperate, meso-oligotrophic freshwater lake. RESULTS: AAP bacteria displayed a clear seasonal trend with a spring maximum following the bloom of phytoplankton and a secondary maximum in autumn. As the AAP bacteria represent a highly diverse assemblage of species, we followed their seasonal succession using the amplicon sequencing of the pufM marker gene. To enhance the accuracy of the taxonomic assignment, we developed new pufM primers that generate longer amplicons and compiled the currently largest database of pufM genes, comprising 3633 reference sequences spanning all phyla known to contain AAP species. With this novel resource, we demonstrated that the majority of the species appeared during specific phases of the seasonal cycle, with less than 2% of AAP species detected during the whole year. AAP community presented an indigenous freshwater nature characterized by high resilience and heterogenic adaptations to varying conditions of the freshwater environment. CONCLUSIONS: Our findings highlight the substantial contribution of AAP bacteria to the carbon flow and ecological dynamics of lakes and unveil a recurrent and dynamic seasonal succession of the AAP community. By integrating this information with the indicator of primary production (Chlorophyll-a) and existing ecological models, we show that AAP bacteria play a pivotal role in the recycling of dissolved organic matter released during spring phytoplankton bloom. We suggest a potential role of AAP bacteria within the context of the PEG model and their consideration in further ecological models.

The high molecular diversity in Noctiluca scintillans is dominated by intra-genomic variations revealed by single cell high-throughput sequencing of 18S rDNA V4.

The application of high-throughput sequencing (HTS) technologies has revolutionized research on phytoplankton biodiversity by generating an unprecedented amount of molecular data in marine ecosystem surveys. However, high-level of molecular diversity uncovered in HTS-based metabarcoding analyses may lead to overinterpretation of phytoplankton diversity due to excessive intra-genomic variations (IGVs). The aims in this study are to explore the nature of phytoplankton molecular diversity and to test the hypothesis. We carried out single-cell metabarcoding analysis of 18S rDNA V4 sequences obtained in single Noctiluca scintillans cells isolated from various sites in coastal waters of China. Results showed that each single N. scintillans cell harbored a high level of IGVs with about 100 amplicon sequence variants (ASVs). The large numbers of non-dominant ASVs identified in N. scintillans cells, which might correspond to the larger numbers of ASVs annotated as N. scintillans and showed similar temporal dynamics in metabarcoding analyses, could inflate the inter-species diversity or intra-species genetic diversity. In addition, there were large numbers of additional ASVs that were not annotated as N. scintillans. These non-N. scintillans ASVs might represent diverse preys for N. scintillans, consistent with previous reports that N. scintillans may act as chance predator of a broad-spectrum preys. This single-cell study has unambiguously demonstrated that the existence of high levels of IGVs in N. scintillans and most likely many other phytoplankton species, demonstrating that the majority of the molecular diversity revealed in metabarcoding analysis, which were generally interpreted as the sum of inter-species diversity and intra-species diversity, actually included high levels of IGVs and should be interpreted with caution.

Plankton community changes during the last 124 000 years in the subarctic Bering Sea derived from sedimentary ancient DNA.

Current global warming results in rising sea-water temperatures, and the loss of sea ice in Arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon sequencing data from sedimentary ancient DNA of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial-interglacial cycles, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by picosized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic picosized chlorophytes and Triparmaceae. By contrast, the glacial period is characterized by microsized phototrophic protists, including sea ice-associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift toward the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community's potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea's function as an effective carbon sink.

Annual recurrence of prokaryotic climax communities in shallow waters of the North Mediterranean.

In temperate coastal environments, wide fluctuations of biotic and abiotic factors drive microbiome dynamics. To link recurrent ecological patterns with planktonic microbial communities, we analysed a monthly-sampled 3-year time series of 16S rRNA amplicon sequencing data, alongside environmental variables, collected at two stations in the northern Adriatic Sea. Time series multivariate analyses allowed us to identify three stable, mature communities (climaxes), whose recurrence was mainly driven by changes in photoperiod and temperature. Mixotrophs (e.g., Ca. Nitrosopumilus, SUP05 clade, and Marine Group II) thrived under oligotrophic, low-light conditions, whereas copiotrophs (e.g., NS4 and NS5 clades) bloomed at higher temperatures and substrate availability. The early spring climax was characterised by a more diverse set of amplicon sequence variants, including copiotrophs associated with phytoplankton-derived organic matter degradation, and photo-auto/heterotrophic organisms (e.g., Synechococcus sp., Roseobacter clade), whose rhythmicity was linked to photoperiod lengthening. Through the identification of recurrent climax assemblages, we begin to delineate a typology of ecosystem based on microbiome composition and functionality, allowing for the intercomparison of microbial assemblages among different biomes, a still underachieved goal in the omics era.

Particle-attached bacteria act as gatekeepers in the decomposition of complex phytoplankton polysaccharides.

BACKGROUND: Marine microalgae (phytoplankton) mediate almost half of the worldwide photosynthetic carbon dioxide fixation and therefore play a pivotal role in global carbon cycling, most prominently during massive phytoplankton blooms. Phytoplankton biomass consists of considerable proportions of polysaccharides, substantial parts of which are rapidly remineralized by heterotrophic bacteria. We analyzed the diversity, activity, and functional potential of such polysaccharide-degrading bacteria in different size fractions during a diverse spring phytoplankton bloom at Helgoland Roads (southern North Sea) at high temporal resolution using microscopic, physicochemical, biodiversity, metagenome, and metaproteome analyses. RESULTS: Prominent active 0.2-3 µm free-living clades comprised Aurantivirga, "Formosa", Cd. Prosiliicoccus, NS4, NS5, Amylibacter, Planktomarina, SAR11 Ia, SAR92, and SAR86, whereas BD1-7, Stappiaceae, Nitrincolaceae, Methylophagaceae, Sulfitobacter, NS9, Polaribacter, Lentimonas, CL500-3, Algibacter, and Glaciecola dominated 3-10 µm and > 10 µm particles. Particle-attached bacteria were more diverse and exhibited more dynamic adaptive shifts over time in terms of taxonomic composition and repertoires of encoded polysaccharide-targeting enzymes. In total, 305 species-level metagenome-assembled genomes were obtained, including 152 particle-attached bacteria, 100 of which were novel for the sampling site with 76 representing new species. Compared to free-living bacteria, they featured on average larger metagenome-assembled genomes with higher proportions of polysaccharide utilization loci. The latter were predicted to target a broader spectrum of polysaccharide substrates, ranging from readily soluble, simple structured storage polysaccharides (e.g., laminarin, α-glucans) to less soluble, complex structural, or secreted polysaccharides (e.g., xylans, cellulose, pectins). In particular, the potential to target poorly soluble or complex polysaccharides was more widespread among abundant and active particle-attached bacteria. CONCLUSIONS: Particle-attached bacteria represented only 1% of all bloom-associated bacteria, yet our data suggest that many abundant active clades played a pivotal gatekeeping role in the solubilization and subsequent degradation of numerous important classes of algal glycans. The high diversity of polysaccharide niches among the most active particle-attached clades therefore is a determining factor for the proportion of algal polysaccharides that can be rapidly remineralized during generally short-lived phytoplankton bloom events. Video Abstract.

Archaeal blooms and busts in an estuarine time series.

Coastal bays, such as Delaware Bay, are highly productive, ecologically important transitions between rivers and the coastal ocean. They offer opportunities to investigate archaeal assemblages across seasons, with the exchange of water masses that occurs with tidal cycles, and in the context of variable organic matter quality. For a year-long estuarine, size-fractionated time series, we used amplicon sequencing, chemical measurements, and qPCR to follow archaeal groups through the seasons. We detected seasonally high abundances of Marine Group II archaea in summer months which correlate with indicators of phytoplankton production, although not phytoplankton biomass. Although previous studies have reported associations between Marine Group II archaea and particles, here they are almost entirely found in very small particles (0.22-0.7 µm), suggesting they are free-living cells. Populations of Nitrososphaeria did not vary with particle size or environmental conditions. Methanogens were significant fractions of archaeal sequences in large particles at low tide during winter months. Contrary to expectations, Nanoarchaeia were found predominantly in the free-living fraction despite the previous observation that they require an association with hosts. These results underscore the utility of time series studies in shallow, tidally mixed estuarine environments that capture variable conditions for understanding the ecology and biogeochemistry of planktic archaea.

Molecular physiology of Antarctic diatom natural assemblages and bloom event reveal insights into strategies contributing to their ecological success.

The continental shelf of the Western Antarctic Peninsula (WAP) is a highly variable system characterized by strong cross-shelf gradients, rapid regional change, and large blooms of phytoplankton, notably diatoms. Rapid environmental changes coincide with shifts in plankton community composition and productivity, food web dynamics, and biogeochemistry. Despite the progress in identifying important environmental factors influencing plankton community composition in the WAP, the molecular basis for their survival in this oceanic region, as well as variations in species abundance, metabolism, and distribution, remains largely unresolved. Across a gradient of physicochemical parameters, we analyzed the metabolic profiles of phytoplankton as assessed through metatranscriptomic sequencing. Distinct phytoplankton communities and metabolisms closely mirrored the strong gradients in oceanographic parameters that existed from coastal to offshore regions. Diatoms were abundant in coastal, southern regions, where colder and fresher waters were conducive to a bloom of the centric diatom, Actinocyclus. Members of this genus invested heavily in growth and energy production; carbohydrate, amino acid, and nucleotide biosynthesis pathways; and coping with oxidative stress, resulting in uniquely expressed metabolic profiles compared to other diatoms. We observed strong molecular evidence for iron limitation in shelf and slope regions of the WAP, where diatoms in these regions employed iron-starvation induced proteins, a geranylgeranyl reductase, aquaporins, and urease, among other strategies, while limiting the use of iron-containing proteins. The metatranscriptomic survey performed here reveals functional differences in diatom communities and provides further insight into the environmental factors influencing the growth of diatoms and their predicted response to changes in ocean conditions. IMPORTANCE: In the Southern Ocean, phytoplankton must cope with harsh environmental conditions such as low light and growth-limiting concentrations of the micronutrient iron. Using metratranscriptomics, we assessed the influence of oceanographic variables on the diversity of the phytoplankton community composition and on the metabolic strategies of diatoms along the Western Antarctic Peninsula, a region undergoing rapid climate change. We found that cross-shelf differences in oceanographic parameters such as temperature and variable nutrient concentrations account for most of the differences in phytoplankton community composition and metabolism. We opportunistically characterized the metabolic underpinnings of a large bloom of the centric diatom Actinocyclus in coastal waters of the WAP. Our results indicate that physicochemical differences from onshore to offshore are stronger than between southern and northern regions of the WAP; however, these trends could change in the future, resulting in poleward shifts in functional differences in diatom communities and phytoplankton blooms.

Assessing the relevance of DNA metabarcoding compared to morphological identification for lake phytoplankton monitoring.

Phytoplankton is a key biological group used to assess the ecological status of lakes. The classical monitoring approach relies on microscopic identification and counting of phytoplankton species, which is time-consuming and requires high taxonomic expertise. High-throughput sequencing, combined with metabarcoding, has recently demonstrated its potential as an alternative approach for plankton surveys. Several studies have confirmed the relevance of the diatom metabarcoding approach to calculate biotic indices based on species ecology. However, phytoplankton communities have not yet benefited from such validation. Here, by comparing the results obtained with the two methods (molecular and microscopic counting), we evaluated the relevance of metabarcoding approach for phytoplankton monitoring by considering different metrics: alpha diversity, taxonomic composition, community structure and a phytoplankton biotic index used to assess the trophic level of lakes. For this purpose, 55 samples were collected in four large alpine lakes (Aiguebelette, Annecy, Bourget, Geneva) during the year 2021. For each sample, a metabarcoding analysis based on two genetic markers (16S and 23S rRNA) was performed, in addition to the microscopic count. Regarding the trophic level of lakes, significant differences were found between index values obtained with the two approaches. The main hypothesis to explain these differences comes from the incompleteness, particularly at the species level, of the barcode reference library for the two genetic markers. It is therefore necessary to complete reference libraries for using such species-based biotic indices with metabarcoding data. Besides this, species richness and diversity were higher in the molecular inventories than in the microscopic ones. Moreover, despite differences in taxonomic composition of the floristic lists obtained by the two approaches, their community structures were similar. These results support the possibility of using metabarcoding for phytoplankton monitoring but in a different way. We suggest exploring alternative approaches to index development, such as a taxonomy-free approach.

Evaluating the use of lake sedimentary DNA in palaeolimnology: A comparison with long-term microscopy-based monitoring of the phytoplankton community.

Palaeolimnological records provide valuable information about how phytoplankton respond to long-term drivers of environmental change. Traditional palaeolimnological tools such as microfossils and pigments are restricted to taxa that leave sub-fossil remains, and a method that can be applied to the wider community is required. Sedimentary DNA (sedDNA), extracted from lake sediment cores, shows promise in palaeolimnology, but validation against data from long-term monitoring of lake water is necessary to enable its development as a reliable record of past phytoplankton communities. To address this need, 18S rRNA gene amplicon sequencing was carried out on lake sediments from a core collected from Esthwaite Water (English Lake District) spanning ~105 years. This sedDNA record was compared with concurrent long-term microscopy-based monitoring of phytoplankton in the surface water. Broadly comparable trends were observed between the datasets, with respect to the diversity and relative abundance and occurrence of chlorophytes, dinoflagellates, ochrophytes and bacillariophytes. Up to 20% of genera were successfully captured using both methods, and sedDNA revealed a previously undetected community of phytoplankton. These results suggest that sedDNA can be used as an effective record of past phytoplankton communities, at least over timescales of <100 years. However, a substantial proportion of genera identified by microscopy were not detected using sedDNA, highlighting the current limitations of the technique that require further development such as reference database coverage. The taphonomic processes which may affect its reliability, such as the extent and rate of deposition and DNA degradation, also require further research.

Environmental DNA metabarcoding reveals the effect of environmental selection on phytoplankton community structure along a subtropical river.

The Dongjiang River, a major tributary of the Pearl River system that supplies water to more than 40 million people in Guangdong Province and neighboring regions of China, harbors rich biodiversity, including many endemic and endangered species. However, human activities such as urbanization, agriculture, and industrialization have posed serious threats to its water quality and biodiversity. To assess the status and drivers of phytoplankton diversity, which is a key indicator of aquatic ecosystem health, this study used Environmental DNA (eDNA) metabarcoding combined with machine learning methods to explore spatial variations in the composition and structure of phytoplankton communities along the Dongjiang River, including its estuary. The results showed that phytoplankton diversity exhibited spatial distribution patterns, with higher community structure similarity and lower network complexity in the upstream than in the downstream regions. Environmental selection was the main mechanism shaping phytoplankton community composition, with natural factors driving the dominance of Pyrrophyta, Ochrophyta, and Cryptophyta in the upstream regions and estuaries. In contrast, the downstream regions was influenced by high concentrations of pollutants, resulting in increased abundance of Cryptophyta. The random forest model identified temperature, dissolved oxygen, chlorophyll a, NO2-, and NH4+ as the main factors influencing the primary phytoplankton communities and could be used to predict changes during wet periods. This study provides valuable insights into the factors influencing phytoplankton diversity and community composition in the Dongjiang River, and demonstrates the application value of eDNA metabarcoding technique in large-scale, long-distance river biodiversity monitoring.

Cell-Penetrating Peptide Delivery of Nucleic Acid Cargo to Emiliania huxleyi, a Calcifying Marine Coccolithophore.

Coccolithophores are a group of unicellular marine phytoplankton that exhibit a prolific capacity for carbon conversion and are critical to ocean biogeochemistry. A fundamental understanding of coccolithophore biomineralization has been limited, in part, by the lack of genetic and molecular tools to investigate the organisms. In particular, it has proven to be difficult to deliver macromolecules across the coccosphere-membrane complex. To overcome this barrier, we employed cell-penetrating peptides (CPP) in the Emiliania huxleyi coccolithophores. We evaluated three established CPPs (TAT, R9, and KFF) and designed a CPP that incorporates a high proline content identified in the protein transduction domain of EhV060, an E. huxleyi virus lectin protein. To measure the delivery performance, we covalently linked CPPs to synthetic peptide nucleic acids (PNA) and attached a fluorescein marker. CPP-PNA-FITC complexes were efficiently delivered across the coccosphere-membrane complex to the cytoplasm of E. huxleyi cells. Characterization of E. huxleyi demonstrates that CPP-PNA are nontoxic and reveals specific effects of CPP-PNA on cell biology and calcification. Direct delivery and characterization of synthetic nucleic acids represent a step forward in synthetic biology to explore coccolithophore biomineralization.

Plankton assessment across the distribution of West African hake and tuna based on eDNA metabarcoding.

The richness of plankton communities determines the fish productivity in the ocean, including important resources that rely on extractive fisheries, such as hakes (genus Merluccius) and tunas (genus Thunnus). Their preys forage on zooplankton, and the latter feed on phytoplankton. Inventories of plankton communities for scientific advice to sustainable fishing are essential in this moment of climate change. Plankton is generally inventoried using conventional methodologies based on large water volumes and visual morphological analyses of samples. In this study, we have employed metabarcoding on environmental DNA (eDNA) samples extracted from small water volumes for plankton inventory from twelve distant sampling stations in the East Atlantic Ocean. Zones rich in hake and tuna prey were detected from eDNA, and multivariate multiple regression analysis was able to predict those zones from diatom-based indices and planktonic diversity based on functional groups. Salinity was negatively correlated with the proportion of diatoms in phytoplankton, highlighting expected impacts of current global change on marine plankton communities. The results emphasise the importance of the plankton richness for fish productivity and support the utility of environmental DNA as a tool to monitor plankton composition changes.

Mock community experiments can inform on the reliability of eDNA metabarcoding data: a case study on marine phytoplankton.

Environmental DNA metabarcoding is increasingly implemented in biodiversity monitoring, including phytoplankton studies. Using 21 mock communities composed of seven unicellular diatom and dinoflagellate algae, assembled with different composition and abundance by controlling the number of cells, we tested the accuracy of an eDNA metabarcoding protocol in reconstructing patterns of alpha and beta diversity. This approach allowed us to directly evaluate both qualitative and quantitative metabarcoding estimates. Our results showed non-negligible rates (17-25%) of false negatives (i.e., failure to detect a taxon in a community where it was included), for three taxa. This led to a statistically significant underestimation of metabarcoding-derived alpha diversity (Wilcoxon p = 0.02), with the detected species richness being lower than expected (based on cell numbers) in 8/21 mock communities. Considering beta diversity, the correlation between metabarcoding-derived and expected community dissimilarities was significant but not strong (R2 = 0.41), indicating suboptimal accuracy of metabarcoding results. Average biovolume and rDNA gene copy number were estimated for the seven taxa, highlighting a potential, though not exhaustive, role of the latter in explaining the recorded biases. Our findings highlight the importance of mock communities for assessing the reliability of phytoplankton eDNA metabarcoding studies and identifying their limitations.

Phytoplanktonic species in the haloalkaline Lake Dziani Dzaha select their archaeal microbiome.

Microorganisms are key contributors of aquatic biogeochemical cycles but their microscale ecology remains largely unexplored, especially interactions occurring between phytoplankton and microorganisms in the phycosphere, that is the region immediately surrounding phytoplankton cells. The current study aimed to provide evidence of the phycosphere taking advantage of a unique hypersaline, hyperalkaline ecosystem, Lake Dziani Dzaha (Mayotte), where two phytoplanktonic species permanently co-dominate: a cyanobacterium, Arthrospira fusiformis, and a green microalga, Picocystis salinarum. To assay phycospheric microbial diversity from in situ sampling, we set up a flow cytometry cell-sorting methodology for both phytoplanktonic populations, coupled with metabarcoding and comparative microbiome diversity. We focused on archaeal communities as they represent a non-negligible part of the phycospheric diversity, however their role is poorly understood. This work is the first which successfully explores in situ archaeal diversity distribution showing contrasted phycospheric compositions, with P. salinarum phycosphere notably enriched in Woesearchaeales OTUs while A. fusiformis phycosphere was enriched in methanogenic lineages affiliated OTUs such as Methanomicrobiales or Methanofastidiosales. Most archaeal OTUs, including Woesearchaeales considered in literature as symbionts, were either ubiquitous or specific of the free-living microbiome (i.e. present in the 3-0.2 µm fraction). Seminally, several archaeal OTUs were enriched from the free-living microbiome to the phytoplankton phycospheres, suggesting (i) either the inhibition or decrease of other OTUs, or (ii) the selection of specific OTUs resulting from the physical influence of phytoplanktonic species on surrounding Archaea.